Spin relaxation of hot electrons in the T- valley of zinc-blende semiconductors 
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We present a technique to calculate the spin relaxation of hot electrons during their energy 
thermalization in the F-valley of zinc-blende semiconductors. The results of this model match 
recent experimental data and they can be used to check the applicability of spin related boundary 
conditions across a semiconductor/ferromagnet junction in reverse bias conditions (spin injection). 

PACS numbers: 



Spin relaxation has an important role in semiconductor 
spintronics devices [l[ . Recent measurements by Crooker 
et al. have shown that the spin polarization across a re- 
verse bias GaAs/Fe junction drops appreciably already 
in moderate bias- voltage levels [2]. Saikin et al. have 
studied the spin relaxation of injected hot electrons for 
relatively large bias conditions and across a wide GaAs 
depletion region Q. The authors have included zone 
edge phonon interactions in a Monte Carlo modeling and 
simulated the spatial spin distribution of injected hot 
electrons in the T, X and L valleys of GaAs. These 
effects were then verified by the electroluminescence of 
Fe- based spin light emitting diodes Ivchenko et al. 
have studied the polarization relaxation of photo-excited 
hot excitons duringa cascade of longitudinal optical (LO) 
phonon emissions [5[ . They have used a method of invari- 
ants to determine the general form of the density matrix. 
None of these models nor the conventional spin relaxation 
description of quasi-thermal electrons [g, 0] is suitable for 
studying the spin relaxation in moderate bias conditions 
0. In this paper, we fill this gap and describe a trans- 
parent cascade spin relaxation process during the ultra- 
fast energy thermalization in the T-valley of zinc-blende 
semiconductors. The energy thermalization is governed 
by emission of long wavelength LO-phonons Q . The re- 
sults of this model match the experimental data and they 
can be used to evaluate the applicability of the boundary 
conditions in modeling spintronics devices. 

The Dyakonov-Perel mechanism Q is by large the 
dominant spin relaxation process during thermalization 
of hot electrons. However, the assumptions that are used 
in calculating the spin relaxation time of quasi-thermal 
electrons [7| are not valid for hot electrons. The bro- 
ken assumptions are that (I) the density matrix of hot 
electrons is far from equilibrium and often is anisotropic, 
(II) the scattering is not elastic, and (III) the spin pre- 
cession of hot electrons due to the intrinsic magnetic field 
is non negligible even within the momentum relaxation 
time (along certain high symmetry directions). 

The spin behavior of hot electrons should be studied af- 
ter each momentum scattering rather than in the asymp- 
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totic limit (after many scattering events). The reason is 
that hot electrons reach the bottom of the conduction 
band after a few LO-phonon emissions and altogether in 
about a picosecond. We begin by writing the intrinsic 
spin dynamics in the T-valley of a zinc-blende semicon- 
ductor @, 



^ = nxs, 

alt 



n 3 -(k) = ah 2 (2m 3 sc E g )- 1 / 2 k j (k 2 e -k 2 m ). 
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The intrinsic Larmor frequency vector, H, has a cubic de- 
pendence in the wavevector components |8j . The {j, £, m} 
subscripts denote any cyclic permutation of {x, y, z}. E g 
is the energy-gap, m sc is the electron's effective mass, 
and a is a dimensionless parameter that relates to the 
strength of the spin-orbit coupling in the semiconductor. 

We assume that the initial energy of a hot electron is 
well above the T point but still below the energy of the 
other valleys in the conduction band. In this case the 
energy thermalization is governed by emission of long 
wavelength LO-phonons via the Frohlich interaction 0. 
After each of the ultra-fast momentum scattering events 
the direction and magnitude of Jl are randomized. The 
probability distribution that an electron with a wavevec- 
tor k will be scattered into wavevector k' is, 



Z(k^k') = 



B{k)8{e{k') + e LO - e{k)) 



k' 2 



2k'kcos(9 k ', k ) 
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The denominator describes the phonon wavevector de- 
pendence of the Frohlich interaction, |<k'|7^F|k>| 2 cx 
|k' — k|~ 2 . We have neglected the weak dispersion re- 
lation between the energy and wavevector in the LO- 
phonon branch, ££o(k' — k) w £lo- B(k) is a normalizing 
factor and it is proportional to the momentum scattering 
time from state k, 
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1 + Vj - e LO /e(k) 
T LO {k) ~ B{k) ~ 87r 2 e(fc) y/l-e LO /e(k) 
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The momentum scattering events are represented by a 
homogeneous Poisson process where the spin distribution 
of a hot electron immediately before its n th scattering 
event is denoted by, 



L(k) 
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where S„_i(k, t) is the solution of Eq. (fTJ) with an initial 
condition S„_ 1 (k, 0) which denotes the spin distribution 
of an injected electron immediately after its (n — l) th 
scattering event. Using the probability distribution in 
Eq. ([3]), the spin distribution immediately after the n th 
scattering event is denoted by, 



S after 1st scattering 



S y after 1st scattering 



S n (k',0) 
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d 3 fcS„_i(k)/(k^k'). 
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This cascade process is complete by defining the initial 
spin distribution of the hot electrons. Due to the impor- 
tance of spin relaxation in spintronics devices, we con- 
sider spin injection from a ferromagnetic contact into a 
semiconductor. Thus, the initial spin distribution is gov- 
erned by the difference of the spin dependent transmis- 
sion coefficients, 



S (k,t = 0)oc|* T (k, V)| 2 -|^(k,VO| 
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where V denotes the reverse voltage bias across the semi- 
conductor/ferromagnet junction. We use a simple effec- 
tive mass model (e(fc) = h 2 k 2 /2m sc ) to calculate the 
spin dependent transmission coefficients due to tunneling 
across a thin triangular-like Schottky barrier jiq |. The 
results of the following discussion are not qualitatively 
affected by the choice of initial distribution. 

The temperature dependence of this cascade process 
is manifested only by a slight modification of the intrin- 
sic Larmor frequency vector (via a change in band-gap 
energy). However, one should recall that LO-phonon ab- 
sorptions may also occur during the thermalization. The 
probability of an absorption process is smaller than the 
probability of an emission process by exp{— SLo/krjT}. 
Thus, in typical room temperature III-V semiconductors 
the absorption processes may somewhat enhance the spin 
relaxation of hot-electrons. In this paper, we ignore ab- 
sorption processes by limiting our discussion to low tem- 
peratures. One can construct, however, a more general 
probabilistic model that incorporates both emission and 
absorption processes in the iterative procedure. 
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FIG. 1: (a) Momentum scattering time in GaAs versus energy 
(in elo =36 meV units), (b) A typical normalized initial spin 
distribution of injected electrons from the Fermi level of Fe 
into GaAs versus the polar angle (with respect to the interface 
normal). The injected spin direction is aligned with the spin 
direction of majority electrons in Fe. 
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S after 3rd scattering 
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S y after 3rd scattering 




S after 5th scattering 



(d) 



S y after 5th scattering 
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FIG. 2: (color online) (a) Spin distribution over electrons' mo- 
mentum polar and azimuth angles after the 1st, 3rd and 5th 
successive LO-phonon emissions. The energy of the injected 
electrons is 0.2 eV above the conduction band F point and 
their spin distribution is normalized and points in the +z di- 
rection (see Fig. [lb). S z keeps its dominant role throughout 
the energy thermalization process whereas S x (not shown) 
and S y are of much lower amplitude. After thermalization 
the electrons have lose ~50% of the net injected S z value (see 
text). 



We use this cascade spin relaxation procedure to simu- 
late the spin relaxation of hot electrons in a 10 16 cm~ 3 n- 
type GaAs. To calculate the intrinsic Larmor frequency 
vector, i7(k) we use E g — 1.519 eV, m sc =0.067, and 
a = 0.07 [11|. To calculate tlo(^) from Eq. we as- 
sume that its minimum is at 130 fs and that £^0=36 meV 
Fig. |TJl shows the energy dependence of the mo- 
mentum scattering time. Fig.[T|D shows a typical normal- 
ized initial spin distribution of injected electrons from 
the Fermi level of Fe versus the angle from the semi- 
conductor/ferromagnet interface normal. The strongest 
spin polarization is from electrons whose motion is along 
the interface normal (cosf? — > 1). This distribution is 
independent of the azimuthal angle (0) that lies in the 
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semiconducotr/ferromagnet interface plane. The injected 
spin direction is aligned with the spin direction of major- 
ity electrons in Fe and we set it as the S z component [l3| . 
Fig. [2] shows the evolution of the spin distribution after 
the first, third and fifth successive LO-phonon emissions 
of hot electrons whose initial energy is 0.2 eV above the 
semiconductor conduction band minima. Assuming that 
these electrons are injected from the iron's Fermi level, 
then at 4K this scenario describes a GaAs/Fe junction 
at -0.2 V reverse bias- voltage (or -0.3 V at 300K). It 
takes five successive LO-phonon emissions for electrons 
to reach the vicinity of the T point. Fig. [5^ shows that 
after the first LO-phonon emission, the main contribu- 
tion to S z is still from electrons whose motion is along 
the interface normal. This can be understood due to the 
forward scattering nature of the Frohlich interaction (see 
the denominator of Eq. ([3])). Thus, the initial distribu- 
tion is not totally randomized in wavevector space after 
the first LO-phonon emission. This is not the case after 
five LO-phonon emissions where we see that the contri- 
bution to the net spin is nearly isotropic in the electron 
momentum. At these injection levels the spin informa- 
tion is largely kept when the electrons reach the bottom 
of the conduction band, and S z is much larger than the 
S x (not shown) and S y (shown) components. 

The important figure of merit is to determine the net 
spin of the injected electrons after their ultra-fast energy 
thcrmalization. Using the effective mass approximation, 
the net spin after n LO-phonon emissions is given by 
k^dk n f dtp J d9 sin (9 5„.i(k„, 0), where i enumerates the 
spin components. This writing assumes that initially the 
electrons' energy is concentrated in a thin energy shell be- 
tween Eq and Eq + deQ. Due to the energy conservation of 
the Frohlich interaction (Eq. ([3])), the ratio between the 
k^dk n and k^dko prefactors simply yields k n /k^. Thus, 
the fraction of the net injected spin that is left after the 
thcrmalization process is denoted by, 



R(e ) 



£0 
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J d(f)fd0 sine S N>z (k N ,0) 
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where Sq is the electrons' injected energy above the semi- 
conductor conduction band T-point and N is the number 
of phonon emissions that are needed to reach the bottom 
of the band. We have considered only the z-direction 
since the injected net spin is aligned along this direction. 



Fig. [3] shows the ratio between the final and initial 
injected net spin. The spin information is largely kept 
when the energy of injected electrons is less than 0.2 eV. 
These results are also in accordance with recent mea- 
sured data by Crooker et al. |2j]. The slower spin decay 
of the experimental results is attributed to the fact that 
the theoretical curve considers only the most energetic 
injected electrons (whose spin relaxation is fastest). In 
the experiment, however, electrons with lower energy also 
contribute to the current (the allowed injection energy 
is anywhere between the conduction band edge in the 
semiconductor and the Fe Fermi level). The application 
of this iterative procedure is limited if the injected elec- 
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FIG. 3: The ratio between the final and injected net spin as 
a function of the initial injected energy. The dots denote the 
measured values by Crooker et al. ci- 



trons experience strong inter-valley scattering processes 
In GaAs this limit is reached when the injection en- 
ergy is £o=0-3 eV [3]. Our model shows that less than 
1/5 of the net injected spin is left during the possible 
n = 8 LO-phonon emissions in this limit. 

In summary, we have presented a transparent tech- 
nique to trace the spin evolution during the relaxation 
of hot electrons in the conduction band T-valley of zinc- 
blende semiconductors. We have shown that if the volt- 
age drop across a reversed biased GaAs/Fe junction is 
moderate then the spin information is largely kept. This 
technique can be used to evaluate the applicability of the 
boundary conditions in modeling spintronics devices. 
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